%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%  EXTRA FUNCTION, SPECIFIED BY USER
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

th = linspace(0,2*pi,5000);
xcirc = cos(th);
ycirc = sin(th);

figure(1);
hold on;
pp1=plot(xcirc,ycirc,'r-');
set(pp1,'linewidth',2);
hold off;

for i=1:NumPhysNodes
  xn=node(i,1);
  yn=node(i,2);
  
  t1=text(xn,yn,[num2str(i)]);
  set(t1,'color',[1 0 0]);
  set(t1,'fontweight','bold');
  set(t1,'fontsize',18);
end

for i=1:NumPhysElems
  xmid = (node(tnode(i,1),1)+node(tnode(i,2),1)+node(tnode(i,3),1))/3;
  ymid = (node(tnode(i,1),2)+node(tnode(i,2),2)+node(tnode(i,3),2))/3;

  t1=text(xmid,ymid,[num2str(i)]);
  set(t1,'color',[0 0 0]);
  set(t1,'fontweight','bold');
  set(t1,'fontsize',18);
end

figure(2)
hold on;
pp1=plot(xcirc,ycirc,'r-');
set(pp1,'linewidth',2);
hold off;

for i=1:SubNumPhysNodes
  xn=sub_node(i,1);
  yn=sub_node(i,2);
  
  t1=text(xn,yn,[num2str(i)]);
  set(t1,'color',[1 0 0]);
  set(t1,'fontweight','bold');
  set(t1,'fontsize',18);
end

for i=1:SubNumPhysElems
  xmid = (sub_node(sub_tnode(i,1),1)+sub_node(sub_tnode(i,2),1)+sub_node(sub_tnode(i,3),1))/3;
  ymid = (sub_node(sub_tnode(i,1),2)+sub_node(sub_tnode(i,2),2)+sub_node(sub_tnode(i,3),2))/3;

  t1=text(xmid,ymid,[num2str(i)]);
  set(t1,'color',[0 0 0]);
  set(t1,'fontweight','bold');
  set(t1,'fontsize',18);
end

figure(1)